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No simplification is possible for the region near the grid, where | e<f>/KT t \ 
may be large. Fortunately, this region does not contribute much to the 
thickness of the cloud (called a sheath), because the potential falls very 
rapidly there. Keeping only the linear terms in Eq. [1.-13], we have 



d 2 <f> nppg 2 ± 

to d7 w* t, • ,4, 



Defining 




[1-15] 



where n stands for nop, we can write the solution of Eq. [1-14] as 

<t> =^ 0 exp(-|x|/A D ) [1-16] 

The quantity A D , called the Debye length, is a measure of the shielding 
distance or thickness of the sheath. 

Note, that as the density is increased, A D decreases, as one would 
expect, since each layer of plasma contains more electrons. Furthermore, 
A D increases with increasing KT t . Without thermal agitation, the charge 
cloud would collapse to an infinitely thin layer. Finally, it is the electron 
temperature which is used in the definition of A D because the electrons, 
being more mobile than the ions, generally do the shielding by moving 
so as to create a surplus or deficit of negative charge. Only in special 
situations is this not true (see Problem 1-5). 

The following are useful forms of Eq. [1-15]: 

A D = 69(T/n) 1/2 m, Tin°K 

[1-17] 

A D = 7430(£77n) 1/2 m, KT xn eV 

We are now in a position to define "quasineutrality If the 
dimensions L of a system are much larger than Ad. then whenever local 
concentrations of charge arise or external potentials are introduced into 
the system, these are shielded out in a distance short compared with L, 
leaving the bulk of the plasma free of large electric potentials or fields. 
Outside of the sheath on the wall or on an obstacle, V*<f> is very small, 
and rii is equal to n«, typically, to better than one part in 10 . It takes 
only a small charge imbalance to give rise to potentials of the order of 
KT/e. The plasma is "quasineutral"; that is, neutral enough so that one 
can take n< — n e =* n, where n is a common density called the plasma 
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We now come to the problem of diffusion in a fully ionized plasma. 
Since the dissipative term contains the difference in velocities v< - v«, 
it is simpler to work with a linear combination of the ion and electron 
equations such that v> - v e is thfc unknown rather than v* or v e separately. 
Up to now, we have regarded a plasma as composed of two interpenetrat- 
ing fluids. The linear combination we are going to choose will describe 
the plasma as a single fluid, like liquid mercury, with a mass density p 
and an electrical conductivity These are the equations of magnetohy- 
drodynamics (MHD). 

For a quasineutral plasma with singly charged ions, we can define the 
mass density p 9 mass velocity v, and current density j as follows: 

p = mM + n e m n(M + fn) [5-78] 

1 A/vi + mv e 

v * - (riiM v, + n t m v t ) ^ — — [5- 79] 

p M + m 

j m eiriiVi - n e y e ) ^ ne(\i - \ e ) [5-80] 

In the equation of motion, we shall add a term Mng for a gravitational 
force. This term can be used to represent any nonelectromagnetic force 
applied to the plasma. The ion and electron equations can be written 

<?v- 

Mn — = en(E + v> x B) - Vp { + Mng + P u 15-811 
dt 

dv 

mn — = -en (E + v € X B) - Vp € + mn g + Y fi [5-82] 
dt 

For simplicity, we have neglected the viscosity tensor it, as we did earlier. 
This neglect does not incur much error if the Larmor radius is much 
smaller than the scale length over which the various quantities change. 
We have also neglected the (v • V)v terms because the derivation would 
be unnecessarily complicated otherwise. This simplification is more 
difficult to justify. To avoid a lengthy discussion, we shall simply say that 
v is assumed to be so small that this quadratic term is negligible. 
We now add Eqs. [5-81] and [5-82], obtaining 

n — (Mvi + mv<) = en(Vi — v e ) x B - Vp + n(M + m)g [5-83] 
dt 

The electric field has cancelled out, as have the collision terms P ti = ~P« 
We have introduced the notation 

to = to, ^ to. [5-84] 



for the total pressure. With the help of Eqs. [5-78H5-80], Eq. [5-83] can 
be written simply 

p — = jxB-V/>+pg 15-85] 

dt 

This is the single-fluid equation of motion describing the mass flow. The 
electric field does not appear explicitly because the fluid is neutral. The 
three body forces on the right-hand side are exactly what one would 
have expected. 

A less obvious equation is obtained by taking a different linear 
combination of the two-fluid equations. Let us multiply Eq. [5-81] by m 
and Eq. [5-82] by M and subtract the latter from the former. The result 
is 

Mmn — ( Vi - v,) = en (At + m )E + en(m?i + Mv e ) x B - m Vpi 
dt 

+ M Vp t - (M + m)P ei 15-80] 
With the help of Eqs. [5-78], [5-80], and [5-61], this becomes 

^ l (L) = epE - (M + m)nem - m Vp> + M Vp t 
e dt XnJ 

+ en(mvi+ Mv e ) X B [5-87] 
The last term can be simplified as follows: 

m Vi + My e = Af\i +mv e + M(v e - Vi) + m(vj v e ) 

= -v-(Af-w)— [5-881 
n ne 

Dividing Eq. [5-87] by ep, we now have 

E + v x b - t, j = - f — I- ( I) + (M - m ) j x B + m Vfc - M Vp e ] 
ep L e dt \nV J 

15-89] 

The dfdt term can be neglected in slow motions, where inertial (i.e., 
cyclotron frequency) effects are unimportant. In the limit m/Af 0, Eq. 
[5-89] then becomes 

E + vXB = 7?j + — (j xB-Vfc) [5-90] 

en 



FUSION PLASMA ANALYSIS 



Weston M. Stacey, Jr. 

Callaway Professor of Nuclear Engineering 
Georgia Institute of Technology 




I A WILEY-BNTTERSCIENCE PUBUCATIGN 

I JOHN WILEY SONS New York • Chichester ♦ Brisbane • Toronto 

i 



o 
a 



I 



d 



1 

I 



1 



o 

I 

] 

CO 



X) 

s 



I 

g 



% 11 




to 

.2 a 

a 8 

8* 
| § 
* 1 

I 



d- 
o 

f 
M 

a> 
o 
d 

-d 4> 

d TJ 
< d 



i 
• 



d rt 
8 f | 
1 1 1 

1 6 * 

S 3 $ . 



* 



s 



# 

DO 



s 

'G 

I 



3 

9 I 



oo 

o 



INTRODUCTION 



7 



plasma with a charged particle number density of 10 20 m -3 and suppose 
that the electron number density (n e ) in a spherical volume of 10~ 3 m 
radius (r) were to differ by one per cent from the positive ion number 
density (n»). Denoting the ion charge by e and the electron charge by — e, 
the total net charge (q) inside the sphere would be 

4 

q = gflT 3 (rci - n e )e (2.1) 
and the electric potential (<j>) at the surface of the sphere would be 

where e 0 is the permittivity of free space. Plugging numerical values into 
(2.2) yields <f> = 6 x 10 3 volts. Recalling that 1 eV = 1.602 x 10 -19 joule, 
we find that kT - 1 eV when T = 11, 600 K, where k is Boltzmann's 
constant (1.380 x 10 -23 joule/ K). Therefore, a plasma temperature of 
several millions of degrees Kelvin would be required to balance the electric 
potential energy with the average thermal particle energy. 

Departures from macroscopic electrical neutrality can naturally occur 
only over distances in which a balance is obtained between the thermal 
particle energy, which tends to disturb the electrical neutrality, and the 
electrostatic potential energy resulting from any charge separation, which 
tends to restore the electrical neutrality. This distance is of the order of a 
characteristic length parameter of the plasma, called the Debye length. In 
the absence of external forces, the plasma cannot support departures from 
macroscopic neutrality over larger distances than this, since the charged 
particles are able to move freely to neutralize any regions of excess space 
charge in response to the large Coulomb forces that appear. 



2.2 Debye Shielding 

The Debye length is an important physical parameter for the descrip- 
tion of a plasma. It provides a measure of the distance over which the 
influence of the electric field of an individual charged particle (or of a sur- 
face at some nonzero potential) is felt by the other charged particles inside 
the plasma. The charged particles arrange themselves in such a way as to 
effectively shield any electrostatic fields within a distance of the order of 
the Debye length. This shielding of electrostatic fields is a consequence of 
the collective effects of the plasma particles. A calculation of the shielding 
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than W y and making use of Eq. 4-30, finally leads to 



p ^ + p( u • V)u + V P=ZE+jX B+ S} - S?u, 
ot 



(4-38) 



where 



(4-39) 



The collisional friction term does not appear in Eq. 4-38 because 2 o R o =0 
from momentum conservation in elastic collisions. 

To obtain a one-fluid equation for the current density, we multiply Eq. 4-7 
by e g /m a and sum over species. Using the velocity decomposition of Eq. 4-37, 
defining the pressure tensor in terms of ^instead of W, and making use of the 
fact that m e /m„«: 1 for o^e leads to 



m. 



n € e] 



*L + v . ( uj+ju - 2 t/u ) + — V • P e 



n e e e 



= \E+{ 1 ^-)uXB + —jXB 



-*j+— 2 Sl 



(4-40) 



where 



n e e\ 0 



^e a JvC a d 3 v, 



a m ° 



(4^2) 



The quantity -q defined by Eq. 4-41 is the plasma resistivity. The subscript e 
refers to electrons. 

Equations 4-38 and 4-40 can be simplified further when the length (L) 
and time (t) scales of interest are large compared to certain characteristic 
parameters of the plasma. Since the potential of a charge is Debye-shielded, 
the plasma is effectively neutral on a length scale L^X Dt so that the 
quasineutrality approximation, 2-»0, can be used. 

In many situations of practical interest, it can be shown that the distribu- 
tion function is an isotropic MaxWellian plus a small correction on the order 
of 8=r Lj /L. In such cases, since the RHS of Eq, 4-4 vanishes for an isotropic 
distribution, v & and u are of first order in 8. Multiplying Eq. 4-38 by L/p, the 
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Macroscopic Plasma Description 



C. MAGNETOHYDRODYNAMIC MODEL 

A specialization of the one-fluid model that is frequently used, together with 
Maxwell's equations, for the analysis of plasmas is known as the magnetohy- 
drodynamic (MHD) model. The MHD model consists of 

+ V • pu= 0, continuity ; (4-49) 

at 

i\ — • 

p-£- +Vp=fxB, momentum ; (4-50) 

tj/= E+uXB, Ohm's law ; (4-5 1 ) 

Vx£ = _H| (4-52) 
dt }, Maxwell's equations; 



(4-53) 



and a suitable "equation of state" instead of the energy balance equation 
such as 

=0, adiabatie; (4-54) 



or 



or 



V • u = 0, incompressible ; (4-55) 



isothermal. (4-56) 



From our previous derivation of the one-fluid equations, we can conclude 
that certain assumptions are implicit in ;the MHD model. The use of a scalar 
pressure assumes an isotropic distribution to lowest order in S=r Li /L. A 
source-free problem is obviously assumed. Neglect of 2 implies that the 
model is valid for spatial scales L^>X D , and neglect of certain other terms 
indicates that L^>r Li is a requirement for validity. The use of a one-fluid 
model implies that the electrons and ions respond together. This implies the 
model is valid on time scales that are long compared to the inverse of the 
lowest of the plasma characteristic frequencies discussed in Chapters 1 and 2, 
which is the ion gyrofrequency; that is, T^flf 1 . In fact, the MHD ordering is 

1»-1-»S 2 . (4-57) 
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FUNDAMENTALS OF PLASMA PHYSICS 



fusion research. 

Another important aspect of plasma behavior is the emission of 
radiation. The main interest in plasma radiation lies in the fact that 
it can be used to infer plasma properties. The mechanisms that cause 
plasmas to emit or absorb radiation can be grouped into two categories: 
radiation from emitting atoms or molecules, and radiation from acceler- 
ated charges. At the same time that ionization is produced in a plasma, 
the opposite process, recombination of the ions and electrons to form neu- 
tral particles, is normally also occurring. As a result of the recombination 
process, radiation is often emitted as the excited particles formed dur- 
ing recombination decay to the ground state. This radiation constitutes 
the line spectra of plasmas. On the other hand, any accelerated charged 
particle emits radiation. The radiation emitted whenever a charged par- 
ticle is decelerated by making some kind of collisional interaction is called 
bremsstrahlung. If the charged particle remains unbound, both before and 
after the encounter, the process is called free-free bremsstrahlung. Radia- 
tion of any wavelength can be emitted or absorbed in bremsstrahlung. 
If the originally unbound charged particle is captured by another parti- 
cle, as it emits the radiation, the process is called free-bound radiation. 
Cyclotron radiation^ which occurs in magnetized plasmas, is due to the 
magnetic centripetal acceleration of the charged particles as they spiral 
about the magnetic field lines. Blackbody radiation emitted from plasmas 
in thermodynamic equilibrium is important only in astrophysical plasmas, 
in view of the large size needed for a plasma to radiate as a blackbody. 

2. CRITERIA FOR THE DEFINITION OF A PLASMA 
2.1 Macroscopic Neutrality 

In the absence of external disturbances a plasma is macro scopically 
neutral This means that under equilibrium conditions with no external 
forces present, in a volume of the plasma sufficiently large to contain a 
large number of particles and yet sufficiently small compared to the char- 
acteristic lengths for variation of macroscopic parameters such as density 
and temperature, the net resulting electric charge is zero. In the interior 
of the plasma the microscopic space charge fields cancel each other and 
no net space charge exists over a macroscopic region. 

If this macroscopic neutrality was not maintained, the potential en- 
ergy associated with the resulting Coulomb forces could be enormous com- 
pared to the thermal particle kinetic enerev. Consider, for examnle. a 
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